Spin waves in ultrathin ferromagnetic overlayers 
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The influence of a non-magnetic metallic substrate on the spin wave excitations in ultrathin 
ferromagnetic overlayers is investigated for different crystalline orientations. We show that spin 
wave dumping in these systems occur due to the tunneling of holes from the substrate into the 
overlayer, and that the spin wave energies may be considerably affected by the exchange coupling 
mediated by the substrate. 
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I. INTRODUCTION 



s ! ■ Recent advances in materials growth techniques. .and precise control of deposition processes have enabled production 
of multilayers with excellent interfacial quality.EfH Presently, it is possible to grow ultrathin magnetic films on a 
'_ ' , substrate, regulating the film thickness very accurately. Epitaxial films with very well-defined thicknesses, showing 
1 virtually no layer thickness fluctuation over macroscopic distances, have been fabricated.!!! The ability to control film 
thicknesses with such accuracy, together with the freedom of choosing different substrates, with distinct crystalline 
■ orientations, broaden the spectrum of magnetic responses, and make these systems highly attractive for technological 
| applications. 

The magnetic behavior of ultrathin films is strongly affected by spin wave excitations. Hence, the study of spin waves 
*^ ' in these structures is important for understanding their magnetic properties and characteristics. In fact, spin waves 
Ch , in ultrathin films have been extensively studied, both experimentally and theoreticallyQuu At low temperatures, for 
O ' instance, the magnetization reduction is basically controlled by long- wavelength spin waves. Thus, measurements of 
M (T) can provide useful information about spin wave excitations in low-dimensional magnetic strucutures. Generally, 
the energy of a long-wavelength spin wave propagating with wave vector q is given by E = D(q)q 2 , where D(q) is the 
exchange stiffness constant. In a relatively thin magnetic layer, spin waves are excited with wave vectors qt\ parallel 
to_the layer, and for certain crystalline orientations, the spin wave energies may depend much upon the direction of 
C\3 q\\U, due to lattice anisotropics. 

It has also been shownQ that spin wave lifetimes in ultrathin ferromagnetic films can be substantially affected by 
the presence of a non-magnetic substrate. By considering a monolayer of a strong ferromagnet on a surface of a 
non-interacting metallic substrate, Mathon et al.Q showed that spin waves in the overlayer become critically damped. 
Such behavior has been attributed to the decaying of spin waves into electron-hole pairs due to the tunneling of holes 
from the substrate into the magnetic overlayer. The substrate also affects the spin wave energy, and modify the region 
in qu -space where it follows a quadratic behavior. 

Here we pursue these ideas and investigate the influence of a non-magnetic substrate in spin wave excitations in 
ultrathin ferromagnetic overlayers for different crystalline orientations. We consider a monolayer of a strong ferromag- 
net both on (100) and (110) surfaces of a semi-infinite nonmagnetic metallic substrate. By artificially reducing the 
electron hopping between the substrate and the overlayer we explicitly demonstrate that the spin wave dumping 
really due to the tunneling of holes from the substrate into the overlayer, as previously pointed out by Mathon et al 
For the (110) overlayer we have found that the spin wave energy depend upon the gjj direction, but such dependence 
is much less pronounced than previously found for unsuported monolayers. This is due to the enhancement of the 
exchange coupling between the local moments in the overlayer mediated by the substrate, which strongly affects the 
spin wave energies. 

This paper is organized as follows: in section 2 we briefly review the theory we have used to calculate the spin wave 
energies and lifetimes in overlayers. In section 3 we present our results and discussions, and finally, in section 4, we 
draw our main conclusions. 
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II. TRANSVERSE SPIN SUSCEPTIBILITY 

The spin- wave spectrum of itinerant ferromagnets can be obtained from the dynamical transverse spin susceptibility 

/oo 

J -°° 

where — {t); S~(0))) is the time-dependent transverse spin susceptibility in real space, given by the 

two-particle Green's function 

xt( t ) = - ie W S t(t)^i(0)]). (2) 

Here, (S^) is the spin raising (lowering) operator at site i, the brackets denote a commutator, (...) represents the 
thermodinamical average which at zero temperature reduces to the ground-state expectation value, and Q(t) is the 
usual step function given by 



e(t) 



if * < 

1 if t > 1 ' 



The exact calculation of xtj (0 involves the solution of an infinite set of coupled equations which in general is not 
possible to find. However, the random phase approximation (RPA) provides an useful decoupling scheme which allows 
one to solve such equations and to obtain the spin wave energies rather accurately. To find an expression for x« _ (0 



within the RPA, we follow references and define the spin operators 

c+ t 

5 y = a W » ( 3 ) 
where a\ a (ai a ) creates (destroys) an electron with spin a at site i. Using S^j, one may define a generalized susceptibility 

Xi+-H(*) = «Stf(*);^(o)». (4) 

Clearly, the transverse spin susceptibility we wish to calculate is xtf CO — xtijj (0- I* 1 order to calculate it, we consider 
that the electronic structure of the system is described by a simple one-band Hubbard Hamiltonian 

H = ^ tij a ia a 3(r + X! + U i n i\ n i{i ( 5 ) 

ija ia i 

where tij is the hopping integral between sites i and j (tu — 0), ej is an atomic energy level, Ui represents the effective 
Coulomb interaction between two electrons on the same site i, and = a^Cbia is the corresponding electronic 
occupation number. 

In this case, xtjki CO obeys the following equation of motion: 

— (t) 

n 

+Uj(([a^a n a^a n ] (t);5«(0))) + [a^a^-J (t); S^(0))) (6) 

The appearance of higher order Green functions leads to an infinite chain of coupled equations for xtjki CO that can 
be decoupled using the RPA which consists of replacing 

a h a jv ,a U a it' ~ ( a L a ^') a U a n' - ( a L a it') a U a j°' + ( a U a ie) a L a ^' - { a ki a w) a \a a w ■ ( 7 ) 

Here, the expectation values are evaluated in the Hartree-Fock ground state, where the dynamics of the two spin 
projections are treated independently. As a result, the average of products of two operators associated with opposite 
spins vanishes. Therefore, one obtains 



3 



ih 



Xl3 ^ - = 5{t){a\ ] ai ] 8 jk - a^aj^u) + ^ (5 im t jn - S jn t mi ) x^nkii 1 ) + 

m,n 

+ ^nC 7 ™ (8im(4i a tt) ~ 5 im(aJ T ai|)J • ( 8 ) 

The Fourier transform of equation || then reads 

huJ xtj U {u) = D m + ^ [KijmnXmnkliu) + 4*mnXmnlW ( w ) + J ijmnXtinkl( u )] > ( 9 ) 
mn 

where £), K, J' and J are four-indices matrices defined by 

Dijki = (a^aifSjk - al^kiSu) 
Kijki — Siktji — Sjitki 

Ji 3 ki = SkiUk \ 8ik(a>ii a jl) ~ ^jk(a\^an)^) 

Jm^^idUjin^-Uiimi)) • (10) 



The matrix elements of the product of two of such matrices is given by ( AB ) = J2 mn Aij mn B mn ki- Thus, we may 
rewrite equation |9| in matrix form as 



hcux + ~(u) = D+{K + J + j')x + ~(u)- (11) 
This equation may be also rewritten as 

- x HF (") + x HF {u)Px + -{u) , (12) 

where P = D^ 1 J, and X {u) represents the susceptibility x^ ( w ) calculated within one-electron theory, i.e., in the 
Hartree-Fock approximation. x HF ( U} ) satisfies the following equation: 

fiu$ HF (w) = D+(K + j')x HF (w). (13) 

Therefore, by using equations ^ we find 

XijfciM = Xijkl(u) - xfjmm^WmXmmkli^) • ( 14 ) 
m 

The dynamic susceptibility Xij~W) = Xujj(t) is then given by: 

Xtfi") = X?j F (u) ~'£xg?(f»)U m x£(u,) . (15) 

m 

In the HF approximation, f- and J,-spin electrons are independent, and an electron with spin a at site i is subjected 
to the HF potential 

e«r = e< + Ui{ni-a) . (16) 
It follows that X§ F ( W ) can be expressed in terms of the HF one-electron propagators as 

x£ F H = -<aj T (*)a iT >4(t) + {at iail (t))Gl*(t) , (17) 
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where Gfj(t) is the time-dependent one-particle retarded Green function for electrons with spin a, connecting sites i 
and j, defined by 



G' 



,(t)=~0{t)({a iir (t),al}), (18) 

where the braces denotes an anticommutator. The correlation functions in equation ^ can be written in terms of the 
one-particle Green functions as 



{a] a a la {t)) = ~- j duf{u) 9fmGf?-(w) e^ 4 , (19) 



where 



Sm G% («) = | [Gy («) - GJ* M] . (20) 



After Fourier transforming equation [l?] and using equation [T^ one obtains 

1 



4*» = -± jdu'W) ZmGj^G^' + ^ + ZmG^G^iuj'-u;) , (21) 



where GJ~ (u>) is the advanced one-particle Green function. 

Since we are interested in multilayer structures having translational symmetry parallel to the layers, it is convenient 
to work with a mixed representation by choosing our basis as Bloch sums in a single atomic plane I defined by 

II jee 

Here tp(Rj) denotes an atomic orbital centered at site Rj € £, gjj is a wave- vector parallel to the layers, and N» is the 
number of sites in this plane. Owing to the fact that gjj is a good quantum number, the Hartree-Fock susceptibility in 
such representation is given by 



|,w) = -iy dco'f(u;')^^[^mGl, e (k ll ,u J , )G l u ,(q ll +k h Lo' + cj)+ 



3m G\ t , (fc|, , J) G\7 t (fe|| - q\\ , J - w)] . (23) 
Similarly to equation pi the in-plane dynamic susceptibility satisfies the following equation: 

xtr(q\\,u) =x?e F (q\\,u) - ^Xt£{?bu)U m x+^(q\\,w) . (24) 

m 

It is noteworthy that equation ^ couples the susceptibility matrix elements involving atomic planes with U =/= 
only. Thus, for a magnetic film of finite thickness on a non-interacting substrate, the set of equations ^ can be solved 
in matrix form as 

X + -{q h uj)=[l + X HF {q h o J )UY 1 X HF {q h uj), (25) 

where [7 + x HF U] is a matrix in plane indices having finite dimension equal to the number of atomic planes of the 
magnetic film. 
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a 


1.0 


0.8 


0.6 


0.3 


0.1 


0.0 


m s 


0.275 


0.292 


0.304 


0.316 


0.320 


0.320 



TABLE I: Magnetic moments of the surface layer (in units of the Bohr magneton) calculated for different values of the hopping 
t± — at between the substrate and the overlayer. Except for a — all values of m 3 have been determined self-consistently. 



III. SPIN WAVES IN SOME OVERLAYERS 

We consider a monolayer of a metallic ferromagnet on a surface of a non-magnetic semi-infinite metallic substrate. 
We examine overlayers placed on (100) and (110) surfaces of a simple cubic lattice. The electronic structure of the 
system is described by the Hamiltonian given by equation |^. We take into account hopping between nearest-neighbor 
sites only, and assume that it is the same (t{j = t) both in the substrate and in the overlayer. We set the atomic 
energy levels ti and the effective on-site Coulomb interactions Ui both equal to zero in the substrate, and fix the 
Fermi energy at Ep = 0.15, which gives a bulk substrate occupancy of n = 1.06 electrons/atom. Here all energies are 
measured in units of the nearest-neighbor hopping t. The HF ferromagnetic ground state of the system is calculated 
self consistently. We choose Ui = 12 in the surface layer, and determine its atomic energy level e s so that the overlayer 
has electronic occupancy of n s — 1.68 electrons/atom (appropriate for Co in this single-band model). A relatively 
large value of U in the overlayer was chosen to guarantee a stable ferromagnetic HF ground state with a number of 
holes in the majority-spin band of the overlayer much smaller than in the minority one. The spin wave spectrum is 
obtained by calculating the surface transverse spin susceptibility x H (<1\\> w )j using equations ^ and ^3| 

First we study a ferromagnetic monolayer on a (100) surface. Figure 1 shows Imx H (gii, ui) calculated as a function 
of energy E = Hlu for several values of qu along the [100] direction in the surface plane. The lifetime of a spin wave with 
wave vector qu is inversely proportional to the width of the peak of Im^ -1 {q\\ , ui) . Such width may be influenced by 
the small imaginary part r\ usually added to the energy in numerical calculations of x" 1 (q, ui). Our calculations for the 
[100] direction were all made with the value of r\ — 1 x 10~ 2 , and for the [110] direction we have used r\ = 5 x 10~ 3 for 
numerical convergency reasons. It is clear from figure 1 that the spin waves in the overlayer become strongly dumped 
for increasing values of qn . The spin wave energies E(q\\ ) were obtained from the position of the peak in Imx H (q«,u>). 
The inset in figure 1 shows that E varies quadratically with qu over a wide range of values in the first Brillouin zone. 

The explanation for the spin wave dumping in the overlayer, given in references |7j], is based on the tunneling of 
holes from the non-magnetic substrate into the majority-spin band of the overlayer. Without a substrate, the free- 
standing ferromagnetic monolayer would have a well-defined Stoner gap, and the spin waves infinitely long lifetimes. 
The tunneling of holes from the substrate into the overlayer destroys this well-defined Stoner gap allowing the spin 
waves to decay into electron-hole pairs. To prove that such explanation is correct we have gradually disconnected 
the overlayer from the substrate by artificially reducing the hopping t± = at between the surface layer and the 
substrate. We consider several values of < a < 1, recalculating in each case the HF ground state self consistently. 
The corresponding values of the surface layer magnetic moments are listed in table 1. As expected, the ground state 
magnetic moment slightly increases as a decreases. Figure 2 shows Imx ^(E, qu), calculated as a function of E for 
different values of a, and q^ — 0.083 x 2n/a, where a is the lattice constant. When a decreases so does the tunneling 
of holes from the substrate into the overlayer. Thus, the decaying probability of the spin waves reduces and their 
lifetimes increase. Consequently, the widths of the spin wave peaks become narrower, as evidenced in figure 2. It is 
also noticeable from figure 2 that the spin wave energies become smaller as a decreases. This is partially due to the 
reduction of the exchange coupling between the local moments mediated by the substrate when a decreases. This 
shows that the presence of a non-magnetic substrate may substantially affect the spin wave spectrum of the monolayer. 

We now examine a ferromagnetic monolayer on a (110) surface. Figures 3 and 4 show results of Im^" 1 (q\\,E), 

calculated as a function of energy for several values of q\\, along the [001] (z) and [110] (£) directions in the surface 
plane, respectively. The spin wave energy for the (110) surface is not isotropic in g||-space. The origin for such 
anisotropy lays on the crystalline structure of the (110) overlayer. For a simple cubic lattice it is formed by chains of 
nearest-neighbor sites along the z direction that in the absence of second-neighbor hopping are linked to each other 
via the substrate only. Thus, without a substrate (i.e., for a free-standing (110) monolayer), those chains would be 
uncoupled, and no energy would be required to excite long wavelength spin waves propagating perpendicularly to the 
chainsB The inset in figure 3 shows that the dispersion relation for spin waves propagating along z varies quadratically 
with q z over a wide range of values in the first Brillouin zone. In contrast, the inset of figure 4 shows that for g| 
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perpendicular to the chains, E(q^) deviates from the quadratic behavior for relatively low values of q^. 

By comparing the insets of figures 3 and 4, we note that the energy of a spin wave propagating along [110] in the (110) 
over layer is smaller than when it propagates along z-direction with the same gjj . The difference in energies, however, 
is not as large as previously found for unsupported monolayers .B The reason is the exchange interaction between the 
chains, that is mediated the substrate. Although smaller, it is comparable to the direct exchange interaction between 
nearest-neighbor sites along the chains, leading to a stiffness along [110] which is of the same order of magnitude of 
that along z. 

By reducing the hopping from the overlayer to the substrate, the inter-chain coupling decreases and the energy to 
excite a spin wave propagating along [110] becomes substantially smaller. This is illustrated in figure 5, and it is 
qualitatively in accordance with what has been found in reference ^. 

IV. CONCLUSIONS 

We have investigated the influence of a non-magnetic metallic substrate on the spin excitations in ultrathin ferro- 
magnetic overlayers. Both the spin wave energies and lifetimes have been determined by calculating the transverse 
dynamic spin susceptibility \^ (<Z|j j E) for different surface crystalline orientations. We have found that the spin waves 
in the overlayer is strongly dumped due to the presence of the substrate. By gradually reducing the hopping between 
the substrate and the overlayer we verify that such dumping is caused by the tunneling of holes from the substrate into 
the overlayer, as previously pointed out by Mathon et aL.H For the (110) overlayer we have found that the spin wave 
energy depend upon the direction of the wavevector q~u with which it is excited. Such dependence, however, is much 
less pronounced than previously obtained for unsuported monolayers. We argue that this is due to the enhancement 
of the exchange coupling between the local moments in the overlayer mediated by the substrate that strongly affects 
the spin wave energies. 
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FIG. 1: Spin wave spectrum for the (100) overlayer. The figure shows Imx H {q\\ , w) calculated as a function of energy E — huj 
for several values of <jj| = q x x along the [100]-direction. Solid line is for q x — 0.042, dashed line for q x — 0.063, dot-dashed line 
for q x = 0.083, and long-dashed line for q x — 0.166. All values of q x are in units of 2ir/a where a is the lattice constant. The 
inset show the corresponding spin wave energies obtained from the positions of the peaks of Imx H . 



FIG. 2: Imx h (gy , E), calculated as a function of E = Tlu), for a fixed q\\ = 0.083 along the [100] direction, and different values 
of the hopping t± = at between the (100) surface layer and the substrate. Thin-solid line is for a = 1.0, dotted line for a = 0.8, 
dashed line for a = 0.6, long-dashed line for a = 0.3, dot-dashed line for a = 0.1, and thick-solid line for a — 



FIG. 4: Spin wave spectrum for the (110) overlayer. The figure shows Imx + (q\\, u>) calculated as a function of energy E = Tiw 
for several values of = g 4 £ along the [llO] direction. Solid line is for q^ = 0.1, dashed line for q% — 0.2, and long dashed line 
for q^ = 0.3. All values of q$ are in units of 2n/a where a is the lattice constant. The inset show the corresponding spin wave 
energies obtained from the positions of the peaks of Imx H • 



FIG. 5: Imx~ + (g||, E), calculated as a function of E = hu), for a fixed value of qn = 0.2^, where £ is a unit vector along the 
[llO] -direct ion in the (110) surface plane. The solid and dashed lines correspond to different values of the hopping t± — at 
between the (110) surface layer and the substrate. Solid line is for a — 1.0, and dashed line for a = 0.1. 



FIG. 3: Spin wave spectrum for the (110) overlayer. The figure shows Imx" 1 (<?]|,^) calculated as a function of energy E = Tito 
for several values of q\\ = q z z along the [001]-direction. Solid line is for q z = 0.1, doted line for q z = 0.2, dashed line for q z — 0.3, 
and long dashed line for q z — 0.4. All values of q z are in units of 27r/a where a is the lattice constant. The inset show the 
corresponding spin wave energies obtained from the positions of the peaks of Imx H 
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